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We consider terahertz absorption and gain in a single miniband of semiconductor 
superlattice subject to a Dichromatic electric field in the most general case of commensurate 
frequencies of the probe and pump fields. Using an exact solution of Boltzmann transport 
equation, we show that in the small-signal limit the formulas for absorption always contain 
two distinct terms related to the parametric and incoherent interactions of miniband 
electrons with the alternating pump field. It provides a theoretical background for a control 
of THz gain without switching to the negative differential conductivity state. For pedagogical 
reasons we present derivations of formulas in detail. 



I. INTRODUCTION 



We consider terahertz absorption and gain in a single miniband of semiconductor 
superlattice subject to a bichromatic electric field in the most general case of commensurate 
frequencies of the probe and pump fields. Using an exact solution of Boltzmann transport 
equation, we show that in the small-signal limit the formulas for absorption always contain 
two distinct terms related to the parametric and incoherent interactions of miniband 
electrons with the alternating pump field. It provides a theoretical background for a control 
of THz gain without switching to the negative differential conductivity state. For pedagogical 
reasons we present derivations of formulas in detail. 

There exists a great interest to generation, amplification and detection of coherent 
terahertz (THz) radiation in semiconductor microstructures. M It is stimulated by a rapid 

u n 

progress of THz sciences and technologies ranging from the astronomy to medicine. [2j, [3| 
Currently, the main direction in this active area of research is the development of quantum 
cascade lascrs| l| and amplifiers. [5j However, these devices require relatively low temperature 
of operation. 

In theory, an intraminiband absorption of THz fields in dc biased semiconductor 
superlatticesp, Q,^ (SSLs) can be negative even at room temperature. It was predicted that 
this high-frequency gain arises in conditions of negative differential conductivity (NDC).[s] 
However, a realization of THz amplifier or oscillator, based on the NDC effect in SSL, is 



complicated by a development of space-charge instability 



111 ] resulting in a formation of 



2 



high-field electric domains inside SSL. This NDC-induced electric instability is destructive 
for THz gain. To suppress the electric instability it was suggested a use of new types of 

12] a stack of short superlatticesfl3| or miniband 



SSLs, in particular lateral superlattices. 

nn 

engineering. [14 Il5l | 

Alternatively, a high-frequency gain can be achieved in SSLs driven by strong ac fields 
(Refs. Negative absorption of aprobe field arises if its frequency and the frequency 



of ac pump field are commensurate. 



16, 



170 In recent Letters [ij, [20|, we clarified parametric 



nature of this effect and found that it can exist without switching to NDC in the time- 
average voltage-current characteristic of SSL. This theoretical result allows to expect that 
the undesirable electric instability can be effectively suppressed in the case of parametric 
gain. 

Here we want to provide some mathematical details which have been omitted in the 
letter 0| . Thus the aim of the present contribution is to derive the formulas describing 
absorption and gain of a THz probe field in biased SSL in the presence of a strong THz 
pump field. Within the semiclassical approach we consider the most general case of arbitrary 
commensurate pump oj\ and probe u>2 frequencies and take into account effects of the relative 
phase (p between the fields. We find that in the limit of a weak probe field the absorption 
can be represented as a sum of the phase-dependent parametric and phase-i ndep endent 
incoherent terms. In particular, our analysis provides a proof of the statement [20| that a 
small-signal net gain of a probe field, which is not corrupted by an odinary generation 
of harmonics, can exist in unbiased SSL only at even harmonics and in biased SSLs only 
at half-integer harmonics of the pump electric field. Absorption at fractional harmonics is 
also considered. Additionally, we compare our formula for small-signal absorption at even 
harmonics with the corresponding formula derived by Pavlovich 1^| in the particular case of 
unbiased SSL and ipo = 0. We demonstrate that his formula formally contains an extra term, 
which nevertheless is identically zero. Therefore, with account of this observation, we can 
confirm that the Pavlvich formula is correct. For pedagogical reasons we present derivations 
of the main equations in this paper in detail. 
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II. HIGH-FREQUENCY ABSORPTION OF AN ARBITRARY PROBE 



We suppose that the total electric field E(t) acting on electrons is a sum of the pump 
E p = E + Ei cos^t) (E is the dc bias) and probe E pr = E 2 cos(uj 2 t + V?o) fields. In a real 
device, E pr may be a cavity mode tuned to the desired THz frequency. We consider the most 
general case of commensurate frequencies 

^ = -, (i) 

uo 2 rn 

where n and m are integers and n/m is an irreducible fraction. Dynamics of the electrons 
belonging to a single miniband is well described by the semiclassical based 
on the use of Boltzmann transport equation 



together with the tight-binding dispersion relation 

1 — cos 



e( P ) = j 



pd 



Here A is the miniband width, d is the SSL period, p is the quasimomentum, f eg 



— - — exp ( — cos^ ) is the equilibrium distribution function lllfl , I is the modified 

2nhiQ \2kbT n J 

Bessel function of the argument A/2/csT. The electron velocity, averaged over the 
distribution function f(p,t), can be calculated as 

V(t) = j V(p)f(p,t)dp, (3) 

where the integration is performed over Brillouin zone \p\ < nh/d, V = de(p)/dp = 
Vosm(pd/h) is the electron velocity, Vo = Ad/2h is the maximal electron velocity in the 
miniband. We define the absorption A(u 2 ) of the probe field E pr (t) as 

T 

1 f^, 



A(to 2 ) = (V(t) cos(cu 2 t + (p Q )) t = - j V(t) cos(uj 2 t + <p Q ) dt, (4) 







where the time-averaging is performed over the common period T = 2imjiO\ = 2ixm/uj 2 of 
both fields. Note that the absorption depends on the relative phase <po- Gain at frequecy u 2 
corresponds to A{u 2 ) < 0. Power P absorbed (P > 0) or emitted (P < 0) in the miniband 
at frequency ui 2 is proportional to A(uj 2 )E 2 . 



The Boltzmann equation ((2j) allows an exact solution (see Appendix A). Using this 
solution we found the expression for absorption in the form 29fl 

A(lU 2 ) = y] y ] 4(ft)4(ft)4-jm(A) [Jh+jn-lifc) + Jl 2 +j n +i(/h)] 



X 



£ JiM J iMJh 

h ,l2=— oo j=— oo 

(O + ^i^i + cos(jmp ) + sin(jny? ) 



(5) 



1 + (O + hwi + kuJ2) 2 r 2 
where J n (x) is the odinary Bessel function, O, = edEi/h (i = 0,1,2) and $ = Oj/tjj 
(i = 1, 2) and A is measured in the units of peak velocity V p = V Ii(y)/2I Q (y) (y = A/2k B T) 
corresponding to the peak current in Esaki-Tsu voltage-current characteristic 0, Q, [a]. The 
derivation of Eq. (jSJ) is presented in Appendix A. 



III. HIGH-FREQUENCY ABSORPTION OF A WEAK PROBE 

In what follows we will consider the small-signal limit E% <C E\. In this limit, we need to 
take only certain combinations of indexes of Bessel functions in j5]). First, we consider the 
most important case when the probe field is a harmonic of the pump field. 



A. Absorption at harmonics {0J2 = mwi) 

In the limit of weak probe field (f3 2 <^ 1), we need to take only the following combinations 
of indexes of Bessel functions in ((5|): (l 2 = 0,j = ±1), (7 2 = 0, j = ±2), (Z 2 = ±1, j = T-2), 
(/ 2 = ±1, j = 0). As a result we represent the absorption as a sum of three terms 

A = J harm (<^o) + (j(A coh + A [ncoh ) + 0{(3l). (6) 

Here J harm depends on the pump amplitude and relative phase <^o as 

(fl + Iuji)t cos ipo + sin ip 



l=—oo 



1 + (O + l^l) 2 T 



2 



(Sl Q + lUl)TCOS<po-8m(p 
+ Jl+m{h>l)- 



(7) 



1 + (O + liO^T 2 

Since J harm does not depend on a small amplitude of the probe field this term is leading in 
Eq. (E]). Note that J harm is just the part of current at m-th harmonic which in phase with 
the probe field. Importantly, negative sign of J harm does not automatically mean a negative 
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absorption of the probe field and P < 0, in contrast to the statement in Ref. [27|. Net gain 
can arise only if other terms in the expression ([6j) are negative. 

Next, A coh describes the parametric amplification of the probe field 2(3] due to its coherent 
interaction with the pump field 



.4 



coh 



= E J ^ 

l=— oo 
— Jl+2m{Pl) 



Jl 



l-2m\h>l 



(Pi 



(fl + Iuji)t cos(2<£> ) + sin(2<£> ) 



1 + (fi„ + loJifr 2 
(VLq + Iuj\)t cos(2<£> ) — sin(2v? ) 

1 + Wo + ^l) 2 T 2 

sin ( 2cp ) 



2 £ J / (/3i)Jz +2m (A) I 



Z=— oo 



+ (fi + Zwi + cu 2 ) 2 r 2 



(8) 



and y4 m describes nonparametric effects 



.4 



incoh 



Z=— oo 



(Q + Iuji + cj 2 )t 



(Oq + lui - u; 2 )r 



1 + (Oq + h)! + 0O 2 ) 2 T 2 1 + (O + lu X - UJ 2 ) 2 T 2 



(9) 



It is important to notice that the sign of A mcoh plays an essential role in stabilization of 



space-charge instability in SSL. 



20 



2l| Namely, the necessary condition of electric stability 



is an operation at a positive slope of time-average current-voltage characteristic. As a rule, 
this slope is positive if A mcoh > 0. Therefore it is desirable to work in conditions which 
provide A incoh > 0. 



1. Unbiased superlattice 

In the case of unbiased SSL (E = 0) the harmonic component of the current becomes 

luJiT cos y? + srn fo 



J harm (^ ) = [1 " ("I)™] E MPl)Jl-m(P 



1 + l 2 U0 2 T 2 
l=-oo 1 

For odd m the contribution of J harm is larger than the sum of the coherent and incoherent 
components of absorption. Therefore, the generation of odd harmonics can blur out a weaker 
effect of parametric gain. In contrast, a generation of even harmonic is forbidden due to 
symmetry, J harm (0o) = 0, and thus for even m 

A = ^(A coh + A incoh ) + 0(/? 2 2 ), (10) 
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where 



A coh 9 V- r (R \j in ^ ^iTcos(2y ) +sin(2y? ) 

A = 2 J APl) J l-2mKPl) 1 + 



Z=— oo 

oo 



Z=— oo 



and 

" 2 ; Z. J '^)l + (^ 1 + £J2 )2 r 2- ( 12 ) 
«=— oo 

In our recent letter|20| we showed that the coherent absorption can be negative (A coh < 0) 
at even harmonics due to a parametric resonance caused by Bragg reflections of miniband 
electrons. Moreover, |v4 coh | has a maximum at an optimal phase y? = y^f^. In this case, 
the total absorption A (flOl) is determined by the concurrence of A coh < and A mcoh > 
and can be also negative (net gain) if |v4 coh | > A mcoh . 

We conclude this section noticing that Eqs. (fill and (fT2l ) have been first obtained by 
Pavlovich. [l|| The relations between our equations and formula of Pavlovich are discussed 
in Appendix B. 



B. Absorption at integer half-harmonics (ll>2 = mwi/2) 



In analogy with the preceding case, in the limit /J2 C 1 we take only the certain 
combinations of indexes of Bessel functions in ((5|): (I2 = 0,j = ±1), (I2 = ±1, j = Tl), 
(/ 2 = ±1, j = 0). As a result we represent the absorption as a sum of two terms 



A{u 2 ) = ^(A coh + A mcoh ) + O(pl). 



(13) 



where 



.4' 



coh 



i=— 00 

— Jl+m(Pl 



Jl-m(Pl 



(O + lui)r cos(2ip ) + sm(2(p ) 



1 + (O + /^) 2 r 2 
(Oo + Iuj\)t cos(2y9 ) — sm(2<^ ) 



1 + (O + l^ifr 2 
2 E Jl(fli)Jl+m(Pi) T~T7n — T~7 1 vTl 

^ 1 + (O + IU) X +UJ2)T 

(=— 00 

1 



1 + (O - h) X - U 2 ) 2 T 2 



sin(2<£>o 



(14) 
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and the incoherent absorption A mcoh is again given by Eq. (J9j) . 

Analyzing (fT4l) for E = we see that A coh = in the unbiased SSL and therefore the 
total absorption becomes 

A(u 2 ) = ^A incoh , (15) 

where A incoh is given by (025 . 

Despite impossibility to have parametric gain at half-integer harmonics in unbiased SSL 
(A coh = 0) , it is still feasible to have a negative absorption without electric instability - i. e. 

A coh < Q) ^incoh > o and |^coh| > ^incoh _ ^ thg cage Q f dc biaged gg L 




C. Absorption at fractional frequencies (u; 2 = muji/n with n > 3) 



First of all let us consider the case n = 3. In the limit ft C 1, we need to take only the 
following combinations of indexes of Bessel functions in ([5j): (Z 2 = 0, j = ±1), (Z 2 = ±2, j = 
=Fl), (Z 2 = ±1, j = 0). As a result we can write the total absorption as a sum of two terms 



A(u 2 ) = ^A incoh + ^-A coh + 0((3l), 



2 

1 incoh 



where the incoherent absorption A m is given by Eq. (J9J) and 

coh ^ 7/ fl \f r (a n (^o + ^i)tcos(3<£>o) + sin(3(^ ) 



.4' 



i=— 00 



(Q + /c<Ji)r cos(3<y?o) — sin(3(y9o) 



+ Jl+m{Pl, 



(=—00 



Jl-m(f3 



1 + (fi + lu X ) 2 T 2 

(O + /cui — 2ci; 2 )t cos(3<^ ) — sin(3v?o) 



1 + (fi + Z Wl - 2cu 2 ) 2 r 2 
(Q + Zcux + 2u; 2 )r cos(3y2 ) + sin(3<^ ) 



(16) 



1 + (fi„ + lu x + 2cj 2 ) 2 t 2 
In particular case E Q = the incoherent absorption is (PT2l) and coherent absorption (TT6 
becomes 



A«* = [1 - (-I)™] £ JlVhW-miPl) 



l=—oo 



Iujit cos(3y?o) + sin(3<£>o) (Zo>i — 2c<j 2 )t cos(3<y2o) — sin(3<^o) 



(17) 



1 + / 2 cu 2 r 2 1 + (lux ~ 2u 2 ) 2 r 2 

Hence, if m is the even number A coh = and the total absorption A is essentially determined 
by the nonparametric term A = j3 2 A mcoh /2. 
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Second, we consider the case of n > 3. In the limit of weak probe /?2 -C 1 one can obtain 
for the absorption 

A(u 2 ) = ^A™ oh + 0(^- 1 )- (18) 

Finally, if the frequencies of the probe and pump fields are incommensurate, the expression 
for the absorption A(ui2) consists of only incoherent component A mcoh . 



IV. CONCLUSION 



In summary, within the semiclassical approach we considered ac current in a superlattice 
driven by a bichromatic electric field. We have obtained the formula for absorption j5]) in the 
most general case of commensurate frequencies of the probe and pump fields using an exact 
solution of Boltzmann equation. We showed that in the small-signal limit the absorption 
can always be represented as a sum of two distinct terms describing the phase-depended 
parametric A coh and phase-independent incoherent A mcoh interactions of miniband electrons 
with the high-frequency pump field. The incoherent component A mcoh describes the free 
carrier absorption modified by the pump. The incoherent component A coh has a parametric 
nature. As follows from Eq. (fT8l) in the case n > 2, the total absorption is dominated by 
the incoherent absorption component. Hence it is possible to get the parametric small-signal 
gain only in the cases when the probe field is harmonic or half-harmonic of the pump field. 

Finally, we would like to notice that effects of nonlinear absorption, including the 
parametric gain, are not limited to semiconductor superlattices and THz electric fields. 
These effects should be observable in different physical systems demonstrating dissipative 
transport in a single energy band. 
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Solution of Boltzmann equation 

In this Appendix we first represent the exact solution of the Boltzmann equation (j2j) in 
the case of total electric field E(t) = E + E 1 cos(u; 1 t) + E 2 cos(u 2 t + <po) , and then derive 
the formula (jS). 

The electron distribution functions / and f eq due to a periodicity in the quasimomentum 
permit a representation in the form of Fourier series 

oo oo 

f( P ,t)= ut)e im *, r{p)= f% eim(p > ( 19 ) 

m=— oo m=— oo 

where tp = pd/H. The Fourier coefficients of the f eq are[n| f%* = dl m (y) /2iThI (y) , I m (y) is 
the modified Bessel function of the argument y = A/2k B T. Substituting (fl9l) in (j2j) we get 



9fm 

dt 



VLq + Vti cos(c<j 1 t) + Q 2 cos(u 2 t + <fo) i m \ fm — (20) 



T T 



Solving the homogeneous differential equation corresponding to ( 1201) , we find 

f m = C(t)e imnot e- t/T exp{im[A sin(^) + (3 2 sin(cu 2 t + <p )}}. (21) 
To find the coefficient C(t) we substitute this expression for f m in (120]) and obtain 

C(t) = fm V 4 W 4 W ; 7 n , ; — ; * exp(-2/ 2 ^ ). 

h,h=—oo v 

Here we used the well-known Bessel formula 

oo 

exp(±z/3sin#) = ^ MP) exp(±iW). (22) 



l=— oo 



Therefore the expression for f m (l2Tj) becomes 

fm = f^ exp{[p 1 sin(u 1 t) + p 2 sin(uj 2 t + <p )]im} 

ii a\r t a\ exp([-il 1 uj 1 - il 2 uj 2 ]t) 
x > JJmft 4 mft - — : r exp ^/ 2 ^ 

h,l2=— oo 
Using again Eq. fl22l) we get 

oo oo 

/m = /m S J h( m l3i)Ji 2 (ml32)Jk 1 ( m Pi)Jk2(mp 2 ) 

l\,l2=— oo ki,k2 =— oo 

exp{[— - &a)u;i - z(Z 2 - A; 2 V 2 ]t} 



1 - i(mQ + hui + 1 2 uj 2 )t 



■exp[—i(l 2 ~ h)<Po\- 
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Let us change indexes Zi — k\ — —ui, Z 2 — k 2 = —v 2 . Then 



oo oo 



fm = fm E E J ll( m A)4M2)4+n( m A)4+^( ra A) 

ll,l2=—oo V\,V2=—<X> 

exp[i(z/iwi + v 2 u 2 )t\ ,. \ /00 , 

1 - 2(miZ + Z^i + Z 2 co> 2 )r 

Now we need to calculate the time-depended average electron velocity (cf. Eq. Q) 

2vr 

F(t) = ^ y y sin(^)/(^)^. (24) 
o 

Substituting the distribution function (1231) into (1241) . we obtain 

^ r fry 

V(t) = -V E f™ / sin (^) exp(zmv?) dip = m—^-(fi - (25) 



Using the definition of the absorption (jlj), (1251) . (1231) and denoting 

1 1 + 7(771^0 + Zicui + Z 2 cj 2 )t 



G(m,Zi,Z 2 ) 



1 - 7(777^0 + Zi^i + Z 2 cu 2 )t 1 + (mfi + hcoi + 1 2 uj 2 ) 2 t 2 ' 



we get 



00 00 



exp(iv 2 ip ) 

x { /1, Z 2 ) - (-1)^+^G(-1, Zx, Z 2 )] (exp^i/iWi + z/ 2 ^ 2 + cj 2 )t]) t exp(^ ) 

+ [G(l, Z 1; Z 2 ) - (-1)^+^G(-1, Z 1; Z 2 )] (exp^z/^i + z/ 2 cu 2 - w 2 )*]) t exp(-^ )} + c.c 

T , y- OO OO 

= ! ^ L1 E E JiM)JiMJ h+ MJi* + M 

® l lt l 2 =— OO V1,V2=— OO 

x {G(l,Zi,Z 2 )(exp[z(i/ 1 c<; 1 + u 2 u 2 + ^ 2 )t]) t exp[2V9 (z/ 2 + 1)] 

- G*(l,Zi,Z 2 )(exp[-z(i/ia) 1 + u 2 u 2 - u; 2 )t]) t exp[-2y2 (z/ 2 - 1)]} + c.c. (26) 

Performing the averaging in Eq. (1261 ) over the common period T = 2im/uji = 2irni/u 2 of 
the pump and probe ac fields satisfying (jTJ), we obtain the conditions 

v-iUi + v 2 oj 2 + u 2 = 0, viuji + /y 2 co> 2 - uj 2 = (27) 

for the first and second term in (l26l) . respectively. It follows from ([271) that z/i = — jm 
(j G Z), as well as that v 2 = jn — 1 for the first term and v 2 — jn + 1 for the second term 
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in (1261) . respectively. As a result we have 

iVoh{y) 



A 



2/o(l/) 



J2 E Jh(Pl)Jh{fr)Jh-jm{fr) 

l\ ,12=00 j = — co 

x [^2+^-1(^2)^(1, Zi, Z 2 ) exp^jn) 

- Ji 2+in+ i(/3 2 ) h, h) exp(-i(p jn)) + c.c. 



Taking into account that V p = Voh(y)/2I (y) with y = A/2/ebT 

i - (O + Ziu;i + / 2 co> 2 )r 



(28) 



and 



iG , (l,Zi,Z 2 )exp(^ jn) 



iG , *(l,Zi,Z 2 )exp(-z< / 9ojn) 



1 + (Oo + hut + l 2 u 2 ) 2 r 2 
i + (O + /icui + Z 2 tu 2 )r 



1 + (O + Z lWl + Z 2 ^ 2 ) 2 r 2 
we obtain the final formula for absorption 



cos(jny9 ) +ism(jnipo)} 



[cos(jnipo) -isin^'n^o)], 



A = V v E E MMMJh-jmiPi) [J h+jn -i{f3 2 ) + J h+jn+1 {i3 2 )} 

ll,l2 = — CO j = — OO 

(O + ^1^1 + Z 2 ^ 2 )TCOs(jn<£> ) + sm(jrnp ) 



(29) 



1 + (O + ZiWi + 1 2 uj 2 ) 2 t 2 
Pavlovich formula 

In our notations the Pavlovich formula for absorption at harmonics and integer-half 
harmonics in unbiased SSL (Eq. (5) in Ref. |l^|) can be written as 



l=—oo 
00 



+ (1/ - l) 2 ( Wl r) 2 2 L 1 ; J 



l=—oo 



ll)J\T 



+ 



{y - Z)(wir) 



(30) 



l + ^r) 2 l + (z/-Z) 2 (cj!r) 2 _ 

where = u) 2 ju\. It is easy to see that the first term in Eq. ( 1301) is the incoherent absorption 
for both v = mj2 and u = m (cf. Eq. (fT2l) ). In the case of absorption at half-harmonics 
{v = m/2) that is the only nonzero term in the expression (l30l) . what is consistent with our 
result Eq. (TT5l) . For absorption at harmonics (v = m) the second term in Eq. (1301) is just the 
coherent absorption (c/. Eq. (fTTj) for </> = 0). Finally we consider the last term in f|30l) for 
v = m. Introducing new summation index k = 2m — I we get 

(m-Z)(o;ir) ^ t/om \ ( m ~ k )( UJ i T ) 



J2 JiWi-M)^ 



+ (to — Z) 2 (o>it) 5 



+ (m — /c) 2 (cj 1 t) 2 



£=— oo k=—oo 

Therefore this term equals to zero, as also easy to check numerically. We conclude that for 
4>o = our equations are identical to the Pavlovich formula. 
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